Cosmic string loops: large and small, but not tiny 
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We develop an analytical model to study the production spectrum of loops in the cosmic string 
network. In the scaling regime, we find two different scales corresponding to large (one order below 
horizon) and small (few orders below horizon) loops. The very small (tiny) loops at the gravitational 
back reaction scale are absent, and thus, our model has no ultra-violet divergences. We calculate the 
spectrum of loops and derive analytical expressions for the positions and magnitudes of the small 
and large scale peaks. The small loops are produced by large bursts of similar loops moving with 
very high velocities in the same direction. We describe the shape of large loops, which would usually 
consist of few kinks and few cusps per oscillation cycle. We also argue that the typical size of large 
loops is set by the correlation length, which does not depend on the intercommutation probability 
p, while the interstring distance scales as p 1 ^ 3 . 



I. INTRODUCTION 

Many cosmological models lead to the production of 
cosmic strings as a result of the phase transition [l|, or 
cosmic superstrings at the end of brane inflation HM 0] ■ 
Very early analytical [3, d, 0] and numerical 0, US] in- 
vestigations agreed that the evolution of infinite strings 
exhibits scaling. An unexpected result came from the 
analysis of string loops by means of numerical simula- 
tions, when the two groups [1, [§] independently found 
that the loops do not scale, and are mostly produced at 
the smallest resolution scale of the simulation. 

To tackle this problem an exact simulation for the cos- 
mic string network in flat space-time [10] was developed. 
The precision of the simulation was only set by unavoid- 
able computer arithmetic uncertainty, which is many or- 
ders of magnitude smaller than the scale of any structure 
in the network. Using this code, it was shown that the 
production of small loops is always accompanied by the 
production of somewhat larger loops [11]. Later, it was 
conjectured that the production of small loops is a 
transient phenomena and in the long run only the large 
loops are produced at the scales of the interstring dis- 
tance ~ O.li. 

The two other groups [T3 | have recently performed 
the numerical analysis with higher resolution versions of 
the simulation codes previously used pip and were able 
to resolve the small loops. Reference [lJ| found some ev- 
idence for a scaling distribution of loops, but their size is 
apparently much smaller than in Refs. (ill. [l2| and Ref- 
erence [lj| found a distribution which diverges at small 
sizes. We believe that the dynamical range of all three 
simulations [1, [H, [H[ is not sufficient, to study the true 
scaling behavior. 

In an attempt to settle the issue, in this article we de- 
velop an analytical model of cosmic strings and compare 



it with another recently proposed model of Ref. [H, [Tj]. 
We show that in the final scaling regime both types of 
loops (small and large) are present, and we do not see 
any evidence of the tiny loops at the gravitational back 
reaction scale. 

The paper is organized as follows. In the next section 
we analyze the production of large and medium-sized 
loops. In the third section we study the production of 
small loops in the vicinity of cusps. We show that the 
dominant scale of small loops is only few orders smaller 
than horizon, and the production of loops at the grav- 
itational back reaction scale is strongly suppressed. In 
the forth section we calculate the spectrum of loops in 
three different regimes with intercommutation probabil- 
ity also taken into account. In the fifth section we de- 
rive analytical expressions for the positions and magni- 
tudes of the small and large scale peaks in the power law 
approximation of the correlation functions. In the last 
section we briefly discuss the effects of cross-correlations 
between opposite-moving waves and compare our results 
with those obtained in Ref. 

II. LARGE SCALES 

The evolution of cosmic strings can be captured by 
right (a(cr)) and left (b(er)) moving waves 

x(a,t) = l( a (a-t)+b(a + tj) 

where x(cr, t) describes the shape of a given string. When 
a self-intersection takes place, the following condition 
must be satisfied 

(a(cr + A - t)-a(a - t)) + (b(<r + A + 1) - b(cr + 1)) = 

where A is the size of a loop that could be formed, and 
t is the time of the intersection. 
For convenience, we define 
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a!(x)dx = a(<j + A)— a(<r) 
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and 



III. SMALL SCALES 



B(cr)=/ b'(a;)da; = b((T + A)-b(cr) 

J a 

The condition of self-intersection can be replaced by an 
inequality 

A(ct - t) + B(cr + t) <e, 

which must be satisfied for some a and t in order for 
a loop of a size from A to A + e to have a chance to 
be formed. In the limit of small e, the self-intersection 
takes place with probability approaching to one. If the 
above inequality holds for an arbitrary small e, then we 
say that a loop is formed with probability p, where p is 
determined from the underlying theory. 

In general, the vectors A(cr) and — B(<r) lie within a 
volume V of a sphere with radius A. However, for the 
extreme scales corresponding to very large (Brownian) 
or very small A, the vectors are confined to a somewhat 
smaller sphere or to a spherical shell respectively. In this 
article we model the production of loops as a random 
process of finding A(cr) and — B(c) such that their dif- 
ference is smaller than e. The intuition tells us, that it 
must be easier to find the right counter-partner if the vol- 
ume V is smaller. However, this is not the whole story 
and other factors must be taken into account in order to 
avoid divergences at small scales. 

If the vectors A(ct) and — B(c) span the entire sphere 
of radius A, then we call loops formed at that scale, the 
medium-sized loops and 

4 3 

Vmedium — 3^"^ ' 

Large loops, are the loops formed from the vectors A(cr) 
and — B(er), which are confined to a sphere of radius cx 

£ \ /^, where £ is the correlation length along the string 

and 

4 3 3 

Vlarge = -TrA^. 

Some of the large loops do not inter-commute with in- 
finite strings throughout their evolution, when other 
quickly rejoin back to the network. As we are only in- 
terested in the loops that are not affected by the infinite 
strings, it is reasonable to expect a sharp cut-off at 



ai ~ 7r£, 



(1) 



for the interstring distance of order £ and intercommuta- 
tion probability p = 1. 

The very large loops are not stable due to their Brown- 
ian shape, regardless of p. After few cycles of oscillations 
the loops fragment into smaller loops, which can only 
decay further by gravitational radiation. Therefore, the 
large loops in the final spectrum (after fragmentation) 
are likely to consist of relatively straight segments with 
at most few kinks (~ 4) and cusps (~ 2 per oscillation 
cycle), which is in agreement with numerical simulations 
ofRefs. fulfil. 



For small loops, the volume occupied by vectors A(cr) 
and — B(er) is approximated by a spherical shell of radius 
A and thickness S with 



V small = ^TrXfS, 



where 



^1 — (-A-)rms 
r a+A 



a {x)dx) rm , 
= (a(er + A) - a(<r)) r 



and 



S=y/(A?)-{\\A\\Y 



Thus, nearly all of the small loops are emitted with very 
high velocities by combining nearly straight segments of 
a- and b- waves pointing in opposite directions. In such 
regions the following condition is likely to be satisfied 
satisfied 

a'(a) = -b'(cr), 

which indicates the presence of a cusp. 

Without loss of generality, we consider a- and b- 
segments of the string such that 

a'(0) = -b'(0) 
a(0) = 
b(0) = 

and expand it in powers of a 

a(a) = CTa '(0) + ^a"(0) + ^a"'(0) + ... 
2 o 

b(cr) = ff b'(0) + yb"(0) + ^b"'(0) + ... 



where 



a'(0)-a"(0)=b'(0)-b"(0)=0 



due to the normality conditions ||a'(cr)|| = ||b'(er)|| = 1. 
At the second order in a, a loop is formed when 



a(ai -t) + b((Ti + 1) = a(<7 2 - t) + b(cr 2 + t) 



or 



(<7i + CT2-2*)a"(0) = (ax+a 2 + 2t)b"(0), 

which is satisfied only if <j\ = a 2 = t. Thus, the small 
loops could not be formed for t ^ at the scales at which 
the third derivative of a(<r) or b(er) is negligible. 
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The relevant statistical quantities are given by 
A 2 = (A'(cr)) rms 

a+A 



= (a'(<r+A)-a'(4 



and 



A (<7 )) rms 
a+A 

= (a> + A)-a») rms . 

The production of loops of size A in the vicinity of a 
cusp is allowed only if the following condition is satisfied 

A 3 A 2 
A 3 (A)— ~A 2 (A) — 

There is a divergence at small scales in the total number 
of loops if 

Urn M > 0, 
a^o 3A 2 

but the total energy emitted by a single cusp remains 
always finite. In the absence of divergences the smallest 
size of loops a s could be approximately obtained from 



3A 2 (q s ) 



0.1, 



(2) 



which correspond to the smallest scale at which the third 
term in the expansion of a(cr) remains of the same order 
as the second term. 



IV. SPECTRUM OF LOOPS 

Without loss of generality, in the reminder of the paper 
we calculate the spectrum of loops at a fixed time 4=1, 
assuming that the full scaling was already established. 
This is done to remove a trivial dependence of the scaling 
production function of loops on time /(y) — > /(A). 

In between the small (a s ) and large (ai) cut-off scales 
there are three different regimes in the production spec- 
trum of loops, which correspond to the populations of 
small, medium and large loops. The production function 
of loops could be estimated as 



/(A)A 2 



A 2 (A) 
V(A) 



A 2 



(3) 



as it must be inversely proportional to the volume V, 
where the vectors A(<r) and — B(cr) live, and directly 
proportional to the rate A 2 at which the volume V is 
explored by these vectors. The two factors of A are due 
to the fact that we are interested in the energy emitted 
into logarithmic bins of the loop length. 



From Eq. [3] we obtain 



/ sma „(A)A 2 oc „ A2(A) ., A 2 



medium 



(A)A 2 oc 



* small (A) 

A 2 (A) 



A 2 (A) 
4tt(5(A) 



2 _ 3A 2 (A) 



V, 



medium 



(A) 



A = 



4?rA 



f (A\A* A 2 (A) A2 3A 2 (A) r- 
fi arge (A)A oc —A = 3 vA. 

Vlarge(A) 4^2 

The numerical analysis of Ref. [18], showed little de- 
pendence of the intercommutation probability p on the 
correlation length £, which is a hidden effect of the cross- 
correlations to be discussed in the last section. Deep in 
the scaling regime, when the vectors A'(cr) are correlated 
over large distances, the structures undergo a multiple 
intersections. This leads to a very high effective proba- 
bility of their intercommutation regardless of p. On the 
other hand, the nearby strings inter-commute with very 
low probabilities, which causes the interstring distance to 
vary with p. Let us assume that 

doc £p K , 

where n is yet unknown parameter. 

In the scaling regime the total energy density in long 
strings scales as 



1 



1 



9 ~ d 2 K e V K ' 

but the effective intercommutation probability remains 
exactly one. Therefore, the probability of a single inter- 
commutation must be proportional to the density of long 
strings, which implies that the production functions of 
loops is 



/(A)A 2 ocp 



- 2re A 2 (A) 
V(A)' 



(4) 



where p can be arbitrary small. 

For very low p, the loops of much large sizes are al- 
lowed to form without ever rejoining back to the infinite 
network. Thus, the large cut-off scale, ai, can be sen- 
sitive to the intercommutation probability. A very large 
Brownian loop of size A would sweep an area with radius 

/a" 2 
1" with about ^ strings passing through the area. 

Each of these strings would intersect the Brownian loop 
roughly 5 times, and thus, 



(5) 



It follows that the total energy emitted into large loops 
is given by 

E total oc p~ 2K / A2[ ^> . A 2 dln(A), 

At the same time, if the decay is dominated by produc- 
tion of large loops, then E to tai oc p~ 2K . Therefore, the 



4 



above integral must not depend on p, which implies that 



k = j and 



d cx £p 3 . 

Remarkably, this result is in a very good agreement with 
numerical simulation of Ref. (la ], where k « 0.31. 

The above analysis shows that the cut-off scale of large 
loops, ai, does not depend on p, and is always given by 
Eq. [TJ The cut-off scale of small loops, a s , is determined 
from the statistical properties of strings which was shown 
not to depend on p. Since the production functions of 



loops scale as p 



P 



" 2 / 3 and the cut-off scales are 



independent of p, the over all spectrum of loops must 
also scale as p~ 2 l z . 



V. EXPONENTIAL APPROXIMATION 

In what follows, we assume that the two-point correla- 
tion function C(u) is approximated by 

C(a) = e- a|CT| . 

The first two parameters which describe the statistics of 
segments A(er) and — B(cr) are calculated exactly 



Ai(A) = 




— - (e-oA-i) 



A 2 (A) = y/2 - 2C(A) = VaA 



\ 3 (A) = ^2a 2 (l~e- A ) 

and the third parameter is approximated in the limit of 
small loops 



2aA 3 



Without the cross-correlations between opposite-moving 
waves taking into account, one can make an estimate of 
the production function of loops in the three regimes 



f smaU (A)A 2 cx 

f medium 

(A)A cx 
fi arge (A)A 2 tx 



P 



-2k /g 



-A" 1 
4tt V 2 



4tt 
3p- 2K 

47T 



VHA-5 ( 6 ) 
\/a£"*A 



The energy scale of small loops can be estimated from 
Eq. [2 



0.15 



There are two more relevant scales a mi and a m2 , which 
specify the range of the medium loops. From Eqs. [6] 
a mi = (6a) _1 and a m2 = £. 
An important ratio is given by 



C; 1 f sm aii{A)A 2 d\n{A) 



total 



where 



Etotai = I f sm au{A)A 2 dln(A) + 
2 / me <w(A)A 2 eUn(A) + 

L 

f large (A)A 2 d\n(A), 

2 

which roughly determines the fraction of energy emitted 
into small loops. For £ = 0.1, a = 3 we obtain 

ai ~ 0.3, 
a, - 0.05 



and 



P - 0.5. 



Thus, about 50% of the energy goes into small loops in 
contrast to 90% obtained in Ref. where the power- 
law approximation of the two point correlation was as- 
sumed. Our model with the power-law ansatz also leads 
to immediate divergences at small scales. In a real net- 
work A3 can in principle be divergent when the kinks are 
present. However, our analysis assumes the Taylor ex- 
pansion only around a given cusp. The probability of a 
near-cusp region to contain recent kinks arbitrary close 
to the cusp is infinitely small. 1 We conclude that phys- 
ically relevant A3 cannot be divergent, and C(a) should 
not be modeled by a power-law. 

The spectrum of loops for the above example is shown 
on Fig. [TJ The qualitative shape of the plot agrees very 
well with numerical simulation of Refs. Never- 
theless, it is too early to expect a full quantitative agree- 
ment between our analytic model and different simula- 
tions, due to a couple of reasons. First of all, at present 
time all of the simulations are not entirely in the scaling 
regime, as the production of loops is dominated by the 
remnants of the initial conditions. Secondly, our model, 
as it stands, is ignorant to the cross-correlations between 
opposite-moving waves. In the following section, we will 
make a first attempt to include such effects into our for- 
malism. 



1 The old kinks are completely smoothed out by gravitational ra- 
diation after a fixed (may be very large) time, but the frequency 
of cusps formation decreases with time. Therefore, in the long 
run, the old kinks can be completely ignored for the analysis of 
small loops formation in the vicinity of cusps. 
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Figure 1: Production function of loops. 

VI. DISCUSSION 

In the limit of negligible expansion the left and right 
moving waves stay practically unchanged throughout the 
evolution. Therefore, 

A(0)-B(<7) >e 

must be satisfied for all < a < £, or otherwise a loop 
of size A would have been formed already. Given that, 
we want to find the probability P that the following in- 
equality holds: 

A(0) -B(0) < e. 

At the moment, let us only consider the small loops. 
We expect P to be smaller, than one would get with- 
out cross-correlations taken into account, as it was done 
in the previous sections. On the other hand, if another 
intercommutation is taking place such that 

A(cr) - B(o-) < e 

for some — £; < a < 0, than P must be much larger. The 
over all effect on the production of small loops will be 
the same as if no cross-correlations are present. 

Although, the production of small loops is a very rare 
process as the correlation length grows, the total effect 



is compensated by a large number of loops produced at 
once. It follows that a number of small loops are likely to 
be emitted in bursts, with very high velocities pointing 
in approximately the same direction. These loops are not 
likely to fragment further, since the third derivative of the 
segments at even smaller scales (< a s ) is not sufficiently 
large, and as a result the small loops will decay only 
gravitationally. 

For the large and medium loops the situation is more 
complicated and at present it is not clear how to take the 
effect of cross-correlations into account. In addition, it 
is desired to include the intersections of nearby strings, 
which inter-commute and produce kinks. All of these 
effects complicate the problem dramatically and can in 
principle lead to some unexpected results. At present we 
simply assume that the production of large and medium- 
sized loops is captured entirely by the statistics of long 
strings, in analogy to the small loops. 

To conclude, we summarize the main differences of our 
results from those obtained in Refs. fl6L [l7| . In the 
scaling regime, we do not expect to see very small loops 
at the scale of the gravitational back reaction. Most of 
the energy is emitted into two different populations of 
the string loops: small and large. The scale of large loops 
is set by the correlation length, which does not depend 
on the intercommutation probability. The scale of small 
loops is determined from the statistics of structures on 
infinite strings, and is expected to be only few orders 
lower than the scale of large loops. For illustration we 
used an exponential form of the correlation function in 
contrast to the power-law approximation used in Refs. 
which cannot correctly represent a physically 
relevant statistics of strings at small scales. 
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